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Abstract. We present a detailed spectroscopic investigation, spanning four winters, of the asymptotic giant branch 
(AGB) star IRAS 02091+6333. Zijlstra & Weinberger (2002) found a giant wall of dust around this star and 
modelled this unique phenomenon. However their work suffered from the quality of the optical investigations 
of the central object. Our spectroscopic investigation allowed us to define the spectral type and the interstellar 
foreground extinction more precisely. Accurate multi band photometry was carried out. This provides us with 
the possibility to derive the physical parameters of the system. The measurements presented here suggest a weak 
irregular photometric variability of the target, while there is no evidence of a spectroscopic variability over the 
last four years. 

Key words, stars: AGB and post-AGB - stars: individual: IRAS 02091+6333 = GSC 04041-01743 



1. Introduction 

When Zijlstra & Weinberger (J2J3(J2J discovered a massive 
dust shell around the AGB star IRAS 02091+6333 they 
could use only a single quick look spectrum and the orig- 
inal TYCHO Bt and Vt magnitudes (Perryman et al. 
1997} to estimate the spectral type and the interstellar 
extinction towards the target. The spectroscopic investi- 
gation was lacking comparison spectra obtained with the 
same instrument setup. Furthermore the TYCHO instru- 
ment was close to its sensitivity limits. Thus we obtained 
spectra and photometry of this unique object for several 
years to derive an accurate spectral type and the fore- 
ground extinction. This allowed us to determine more pre- 
cisely the distance of the target for the modelling of the 
dust shell found on Infrared Astronomical Satellite (IRAS) 
images. Zijlstra & Weinberger (2002) outline such shells 
for various types of objects at late stages of their evolution. 
Their focus was especially on a swept up shell with a void 
in the interstellar matter (ISM) around the target. This is 
crucial for both, the "Swiss cheese" like structure of the 
ISM and for the hydrodynamic evolution of the planetary 
nebula (PN) built after this evolutionary stage. 
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We thus investigated the spectroscopic and photomet- 
ric properties of the central star to provide new input for 
the complete model of this unique object and its surround- 
ings. 



2. Data 

The spectroscopic data were obtained with the Innsbruck 
60 cm telescope (Kimeswenger 120011 Bacher et al. I2001fl 
and an OptoMechanics IOC spectrograph. A CompuScope 
Kodak 0400 CCD camera was attached to the spectro- 
graph. In the year 2000 a grating with 600 1/mm was used, 
resulting in a resolution of 1 A/pixel. Later a grating with 
240 1/mm (2.6 A/pixel) was mounted to achieve a better 
S/N ratio. The complete log of the observations is given in 
Table ^ The exposure times varied with the wavelength 
region from 300 to 900 seconds. Each spectrum is a com- 
posite of individual takes covering about 800 A with the 
high resolution grating and 2000 A with the medium res- 
olution mode. The wavelength shifts were selected to ob- 
tain appropriate overlaps of at least 25% of the wavelength 
range. In most of the nights several spectra were obtained 
to improve the S/N and eliminate cosmic ray events. The 
positions along the slit and the starting wavelengths var- 
ied to eliminate possible systematic errors in the setup. 
Flatfield, bias subtraction and wavelength calibration were 
carried out in a standard manner with the help of MIDAS 
routines. For comparison the late type standard stars 
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HD 23475 (M2IIb), HD 39801 (M2I), HD 13325(M3III), 
HD 40239 (M3II), HD 44478 (M3III), HD 42995 (M3III) 
HD 5316 (M4III) and HD 12292 (M5III) were observed 
with the medium resolution setup. The spectra were not 
calibrated absolutely, but only corrected for exposure time 
and airmass (with the atmospheric extinction given by 
Tug lTMl . 

The wide band direct imaging in Innsbruck was ob- 
tained with an AP7p SITe 502e CCD device (Kimeswenger 
et al. l2()()2al Lederle & Kimeswenger 2003). 24 images were 
taken in the nights of 10 th and 11 th of December 2002 with 
B, V, R and Jc filters. The exposure times were 30, 20, 10 
and 10 seconds in B, V, R and Iq respectively. After stan- 
dard basic CCD reduction the source extraction was per- 
formed using SExtractor V2 (Bertin & Arnouts ll996|l . The 
rms of the comparison standards in the field was < CP01 
(<t b = m 003, cry = m 009, a R = m 007 and 07 = m 004). 
Absolute calibration was obtained in both nights using the 
standard stars HR 580 and HR 596. The rms variations of 
the zero-points were less than m 03, resulting in an error 
of the mean zero-point of less than nl 015. 

Table 1. Log of the Innsbruck 60cm observations 



Date 


mode 


band / 


resolution 








wavelength 


sampling 
A /pixel 


off.* 

A 


05.02.2000 


spectra 


6000-9000A 


1.0 


6.2 


15.01.2001 


spectra 


5000-10000A 


2.6 


10.6 


18.01.2001 


spectra 


5000-10000A 


2.6 


10.6 


11.02.2001 


spectra 


5440-10000A 


2.6 


9.9 


15.02.2001 


spectra 


5000-8860A 


2.6 


8.6 


16.02.2001 


spectra 


5440-8860A 


2.6 


8.0 


06.03.2001 


spectra 


5440-8860A 


2.6 


10.7 


08.01.2002 


spectra 


5600-10000A 


2.6 


13.4 


09.01.2002 


spectra 


5000-9000A 


2.6 


14.1 


10.12.2002 


imaging 


BVRIc 






11.12.2002 


imaging 


BVRIc 






17.01.2003 


spectra 


5000-10000A 


2.6 


17.1 



Effective resolution was measured as FWHM of night sky lines. 

The long term monitoring of the target was done with 
the P7 photometer attached to the 1.2m Mercator tele- 
scope located at La Palma, using the 7 filters of the Geneva 
photometric system (Meylan & Hauck [T9811 Golay^994). 
The data was obtained and reduced within the framework 
of monitoring stars with constant airmass as described in 
Burki et al. (1995). The typical errors derived by the rms 
of the standard stars for the nights are below 0™01. The 
results are listed in Table [3 

3. Results 

3.1. Spectroscopy 

First the spectra were investigated with respect to vari- 
ability. We cannot detect spectral variability (except in 
the region of the telluric water vapor bands longwards to 
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Fig. 1. The average spectrum (upper panel) with the most 
dominant features identified and the deviation from the 
mean spectrum of those spectra, using the same instru- 
ment setup, during winter 2001, 2002 and 2003. For clarity 
an offset of 1.0 was added to each spectrum. No spectro- 
scopic variability in the lines or in the continuum trends 
can be found. 



8000A- see Fig. Q). Also individual equivalent widths of 
metal lines of the IR Call triplet and of Mg lead us to 
the same conclusion. Thus a composite spectrum was de- 
rived for the classification. The classification criteria of 
Kirkpatrick et al. I|1991[) were applied in the same way 
as described in Kimeswenger et al. (2003|. This results in 
a spectral class of M4 - M5. The scheme of TiO indices 
(Malyuto et al. H997[) argues for an M4.2III star. Applying 
the extension by Schmidt-Kaler & Oestreicher Ijl998|l we 
obtain an M4III star with M bol w -2™65 (M v * -0^5). 
The relation of the TiO to the NaD and Mgb indices pro- 
pose luminosity class III. 

The most reliable classification is the direct compari- 
son with standards, using the identical instrument setup. 
For this purpose eight comparison stars were taken during 
February 2001 and January 2003 between the target expo- 
sures with the same setup. Especially the TiO features at 
7200-7400 A and 7700-8000 A are extremely sensitive to 
the effective temperature of the stars in that domain. The 
latter was not covered by the spectrum used by Zijlstra 
& Weinberger H2002fl . The classification using such bands 
is completely independent of the correction for the atmo- 



S. Kimeswenger et al.: Optical Spectroscopy of IRAS 02091+6333 



3 



spheric and interstellar extinction. The results are shown 
in Fig. [5] The target superimposes nearly exactly the 
M4III star HD 5316. The overlay gives us an exact solution 
for the interstellar extinction of E B _ V = (F43 ± (TQ1. 

We tried to use the line ratios of the IR Call triplet 
of the comparison stars and those of the target - all taken 
with the same instrument setup. But not even for the com- 
parison stars, having a much better S/N ratio than the 
target, the effective resolution (see Tabled allow to de- 
rive useful separations between luminosity class III and 
II. 

We favor a spectral type of M4III. The luminosity class 
remains somewhat uncertain, but is most likely III. This 
is supported also by the distance (see section 1331 and dis- 
cussion in Zijlstra & Weinberger 2( 102)1 . 

3.2. Photometry and Distance 

We have searched literature and telescope archives for 
photometry, but found in the optical and in the MIR-FIR 
only those mentioned already in Zijlstra & Weinberger 
PTO^I - namely TYCHO in the optical, Midcourse Space 
Experiment (MSX, Egan et al. W) and IRAS in the 
infrared. The results of our photometry, together with re- 
sults from literature are collected in Table |3J For the flux 
to magnitude conversion for IRAS we used the results of 
Wainscoat et al. I|1992|) and for MSX that of Cohen et al. 
pWijl . The calibration of IRAS by Hickman et al. fiffity 
was not taken into account, as the MSX and the IRAS 
12 /jm fluxes gave differences of about 0™40 for all M type 
comparison stars here. For the TYCHO measurements 
we applied the color equations as given in the catalogue: 
V = Vt-0.090x(B t -Vt) and B-V = 0.850x (B T -V T ). 

During the reviewing process of this work, the 2MASS 
NIR data was made public. Although we listed the values 
in Tabled they were not useful, as the target was overex- 
posed even in the short 53 ms preexposure of the survey. 
We thus didn't include these values in our interpretation. 
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Table 3. Summary of photometric observations except 
those listed in Table [21 The TYCHO measurements were 
converted to the standard system (see text). 
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band 


man 


quality 
or error 


ref. 


10 /II 19 9009 
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9T3 


? 
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I 


7 m 3 


? 


[i] 


5.1.1999 


J 


5 m 326# 


(T03 


2MASS 




H 


4™350 # 


(T29 


2MASS 




K s 


4 m 034 # 


(T34 


2MASS 


» 1995 


[8.3] 


3 m 80 


qal 4 


MSX 




[14.6] 


2^9 


qal 1 


MSX 


« 1983 


[12] 


3 m 76 


qal 3 


IRAS 




[25] 


3 m 55 


qal 3 


IRAS 



[1] Zijlstra & Weinberger J21TO21 

* The 2MASS data is flagged in the data base as 

"radial profile fitting of overexposed source" and 

"profile Et very poor". 



The first TYCHO catalogue (Perryman et al. ITM7|> 
and the reanalysis of the same data in TYCHO-2 (H0g 
et al. 2000) give completely different results for the Bt 
band. The Vr band differs considerably more than the 
error given in the catalogue. Thus we assume the error 
to be significantly higher. For further analysis we use the 
TYCHO-2 values. 

The Mercator monitoring (Fig. |2J clearly shows an ir- 
regular variability with an amplitude below 0™1 in V. 
On the other hand the variation between our measure- 
ments, those of the amateur observations reported in 
Zijlstra & Weinberger ^zWzt) and TYCHO-2 in the opti- 
cal give an amplitude of some tenths of a magnitude. The 



(Vr 



5 ) value (15% and 85% of the 228 photomet- 



Fig. 3. The Mercator photometry shows a weak photo- 
metric variation on medium to long timescales. 



ric transits) given in TYCHO-2 is 1™71. But we assume 
the latter is mainly due to the limits of the experiment. 
The MSX and IRAS measurements together with the spec- 
troscopic stability lead us to the conclusion that the object 
seems to be very stable in the effective temperature. That 
fact limits the optical photometric amplitude to those 
found in the Mercator monitoring. Such small variations 
seem to be usual for late type stars around the tip of the 
first giant branch (Ita et al. (2002 ) . 

The colors (V-R) = l m 15 and (V-Ic)o = 2 m 68 in- 
dicate a M4.5III MK class (Drilling & Landolt HMty . 
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Table 2. Results of the Mercator photometry in the Geneva system (Meylan & Hauckjl981). 



MJD 




V 




u 




b 






ba 


Vl 




.9 


weight 


DZLvo 


Do I OU 


10 


445 


15 


647 


12 


027 


lo 


coo 
ooo 


13 


011 


1 1 




11 


. 281 


4 




ooyuo 


10 


452 


15 


540 


12 


039 


lo 


O'iO 


12 


995 


1 1 


zy i 


11 


.281 


4 


C01 OA 


O I *±U'± 


10 


450 


15 


643 


12 


017 


1 O 


010 


13 


007 


1 1 
1 1 


ZO / 


11 


. 286 


3 


Oz 111 


? u*±oy 


10 


380 


15 


416 


11 


965 


1 O 


OOO 


12 


948 


1 1 
1 1 


zoz 


11 


. 214 


1 




D^ooZ 


10 


427 


15 


637 


12 


030 


lo 


DO / 


13 


035 


1 1 


Zol 


11 


. 267 


2 


52299 


37930 


10 


421 


15 


427 


12 


045 


13 


700 


13 


010 


11 


280 


11 


.261 


3 


52299 


42229 


10 


419 


15 


806 


12 


060 


13 


690 


12 


997 


11 


257 


1 1 


. zoo 


O 


52500 


69875 


10 


458 


15 


606 


12 


055 


13 


659 


13 


082 


11 


304 


11 


.297 


2 


52501 


69090 


10 


453 


15 


472 


12 


047 


13 


676 


13 


067 


11 


306 


11 


.286 


2 


52502 


69127 


10 


447 


15 


714 


12 


047 


13 


680 


13 


022 


11 


303 


11 


.278 


2 


52505 


64131 


10 


434 


15 


573 


12 


027 


13 


634 


13 


034 


11 


276 


11 


.254 


3 


52506 


69780 


10 


421 


15 


552 


12 


026 


13 


650 


13 


018 


11 


279 


11 


.258 


3 


52508 


69154 


10 


430 


15 


629 


12 


040 


13 


676 


13 


017 


11 


262 


11 


.264 


2 



Using our photometry and the color equations above we 
receive V T - [8.3] = 5 m 56. Cohen et al. pT07))l are an- 
nouncing 5'. n 21, 5™58 and 6'. n 02 for M3, M4 and M5 re- 
spectively, which is also consistent with the classifications 
found above. The Mercator monitoring was not used for 
the colors, as this system is lacking investigations, abso- 
lute calibrations and color equations for very late type 
stars (Meylan & Hauck lEMl Moro & Munari HMty . 
We derived, again using TYCHO, MSX measurements and 
HIPPARCOS parallaxes, the spectral energy distributions 
(SED) for M3-5 giants from our spectroscopic sample and 
from those one of van Belle et al. ( 1990$. The latter is a 
sample of stars investigated in detail with interferometers 
for the determination of accurate size and luminosity of 
late type stars. Two stars - namely HD 13325 (M3III) and 
HD 5316 (M4III) - overlap between their and our sample. 
One can clearly see in Fig. 0] the broad band mid-IR flux 
hardly varies between M3 and M4 stars. Only M5 stars 
have a significant different SED. We thus fit our target 
between the M4 stars of the sample. This provides us with 
(using all bands) a distance of 1020 pc. The HIPPARCOS 
errors for the bright comparison stars and the uncertain- 
ties from the photometry of the target results in a lower 
limit of 870 pc and an upper one of 1150 pc. The luminos- 
ity derived from the spectroscopy and the V photometry 
give us 920 pc. According to van Belle et al. (1999) for 
K — [12] > 0™6 there is no correlation with the effective 
temperature anymore. Thus this band, lacking here, would 
deliver us no additional information. 

3.3. Astrometry 

The target was suspected to have a relatively high proper 
motion in the TYCHO catalogue, although the positional 
error given there is very high. We thus tried to derive 
a value for the proper motion in the same way as in 
Kimeswenger et al. Ij2002bjl by using the six sky survey 
plates of the region. But as the target is extremely over- 
exposed and lies near the corner of the POSS-E (1952 and 
1993) and of the Quick- VN (1983) plates, we cannot de- 
rive a proper solution. The accuracy is about the same as 




wavelength [|jm] 



Fig. 4. The SEDs of the M3-5 III stars of the sample from 
van Belle et al. (1999) and of the stars used for spec- 
tral comparison here. The distances for calibration of the 
absolute magnitude were taken from HIPPARCOS. The 
photometry was taken from SIMBAD (By), IRC (K, 
Neugebauer & Leighton lT^|l . MSX (4-21/mi) and IRAS 
(12 and 25/xm). The error bars give the tolerance added by 
the HIPPARCOS errors for the worst case star. The target 
fits well within the two M4III stars assuming a distance 
of 1020 pc. 



in the the TYCHO experiment. Kislyuk et al. I|1999fl give 
similar numbers in their online version, using their plate 
from 1983 and the AC2000 catalogue position from 1905 
(Urban et al. I1997|l . But forwarding the errors in both 
catalogues to the resulting proper motion we again end 
up with large error bars (12±7 mas yr -1 ). The direction 
of the motion vector is co- lined with the TYCHO result, 
suggesting a genuine movement which leads to a motion of 
about 60 km s _1 relative to the surrounding stellar field 
after subtracting the galactic rotation for a distance of 
about 1 kpc (only about 1.5 mas yr -1 ). This is, although 
still possible, very high for galactic disk stars. The vector 
is almost opposite to the galactic rotation parallel to the 
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galactic plane. Further investigations about the popula- 
tion membership are needed. 

4. Conclusions 

Table 4. Summary of the results for IRAS 02091+6333 



spectral type 


M4 ±0 .3 HI 


T c ff 


3350 ±50 K 


type of photometric variability 


irr. or SR 


^ V ^Johnson 


10 m 52 ±0 .o3 


interstellar Eb-v 


O m 43 ±0 .oi 


visual amplitude 


< (Tl 


My 


-(T5 


M bo i 


-2 m 65 


distance 


1020±ilg pc 


luminosity 


920±i5o L Q 



The precise spectroscopic classification of the central 
star IRAS 02091+6333 of the unique giant dust shell re- 
ported in Zijlstra & Weinberger |2002|) allowed us to derive 
a more accurate distance of the target. This leads us to a 
somewhat higher mass estimate of 5.2 M and a size of 
1.8 x 10 17 m (5.2 pc) and thus even enlarges the discrep- 
ancy in lifetime as discussed in section 3.3 of Zijlstra & 
Weinberger (2002) and therefore clearly excludes the sce- 
nario of the mass originating mainly from the star itself. 
This result is strengthened by the fact of the lower lumi- 
nosity we derived (40% below the one estimated there). 
The star lies in the evolutionary tracks (Girardi et al. 
2000a) nearby the red giant branch tip and rather far from 
the TP-AGB. The small irregular photometric variability 
fits well to this assumption (Ita et al. I2002|l . So the model 
of the swept up ISM is strongly supported by our mea- 
surements. The fast movement, although uncertain, may 
shorten the timescales for the bubble. This may suggest a 
wind having a higher momentum. One can only speculate 
on the origin of this additional source of energy. There 
is no signature in the spectroscopy for a hot companion 
causing a symbiotic Mira. Also Schmeja & Kimeswenger 
( 2001) have shown that those objects have clear signatures 
in their NIR and MIR bands. 

IRAS 02091+6333 has a distance of D = 10201^ pc, 
a luminosity of L = 920 ± 150 L Q and, according to 
van Belle et al. I|1999f) . an effective temperature of T e g = 
3350 ± 50 K. Both the spectroscopy as well as the photo- 
metric consistency of the colors with the MK type from 
visual to IR bands suggest there is currently no domi- 
nant circumstellar shell material. Girardi et al. (2000b) 
discussed inhomogeneities in the tracks and in the evo- 
lution due to different He-core burning phases. It affects 
slightly the luminosity (and thus the distance) but hardly 
the photometric colors (Girardi et al. 2000a). This adds 
some uncertainties to the values derived for individual 
sources. 
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